A note on rattlers in amorphous packings of binary mixtures of hard spheres 
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It has been recently pointed out in Refs. [1, 2] that 
our theoretical result for the jamming density of a bi- 
nary mixture of hard spheres [3] "apparently" violates 
an upper bound that is obtained by considering the limit 
where the diameter ratio r = Da/Db — > oo. We believe 
that this apparent contradiction is the consequence of a 
misunderstanding, which we try to clarify here. 

We will follow the notations of our paper [3] and denote 
by 7/ the volume fraction of small particles. Then, a very 
simple argument [1,2] leads to the following result in the 
limit r — >• oo: 
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where ipj t i ~ 0.64 is the jamming packing fraction of 
the monodisperse system (r = 1). In the following, we 
consider a binary mixture of hard spheres undergoing 
some dynamics at fixed temperature ksT =1/(5 (the 
precise value of which is irrelevant and fixes the time scale 
of the problem). The two terms in Eq. (1) correspond to 
very different physical situations that might happen in 
the limit r — > oo: 

• For small enough 77, the large particles form a 
jammed structure, while the small ones diffuse 
freely (as in a liquid) in the voids of the large ones. 

• For large enough 77, all particles participate to the 
jammed structure, in which the large particles are 
"diluted" in a "sea" of small particles. Hence, in 
this regime, most of the large particles only touch 
the small ones, while large - large contacts due to 
density fluctuations are rare. It can be easily seen 
that, according to our theory [3], this qualitative 
picture becomes exact in the limit r — > 00 at fixed 
77 or for 77 — > 1 since, in these limits, the average co- 
ordination number of large spheres around a given 
large sphere at jamming, goes to zero. 



These two regimes are separated by 77^ 



0.265 



where the two terms in Eq. (1) are equal, see Fig. 1. 

In the small cage approximation that we use in our 
theory [3, 4], we assume that a particle is allowed to vi- 
brate within a cage made of the particles surrounding it. 
Moreover we assume that the structure of the system is 
frozen, and the spheres (belonging both to the large and 
small components) can't have any kind of diffusive mo- 
tion in the system. Therefore our theory cannot describe 
the region where the small particles diffuse in the voids of 
the large ones; this happens for large enough r and small 




FIG. 1: The bound for the jamming density in Eq. (1) (thick 
solid lines) and the theoretical value from Ref. [3] for r = 10 12 
at a complexity E = 1.7 (dashed line). The solid lines are 
predictions from our theory [3] for different values of r at E = 
1.7, filled symbols are numerical results from simulations [3] 
and open symbols are experimental data from Ref. [5]. Note 
that the large r - small r\ region is out of the domain of validity 
of our theory, indeed the theoretical curves violate the upper 
bound in that region. 



enough 77. For this reason we avoided the large r-small 
77 region in the comparison with experimental data, see 
Fig. 1 (and the caption of Fig.l of Ref. [3]). 

On the other hand, the jamming packing fraction pre- 
dicted by the theory in the infinite size-ratio limit (shown 
in Fig. 1 for r = 10 12 ) is always below the upper bound 
given by Eq. (1) in the region 77 > r\th where the theory 
is expected to hold, see Fig. 1. 

We performed simulations in order to study the num- 
ber of rattlers in different mixtures. We used the code 
developed by Donev and co-workers [7, 8] based on the 
Lubachevsky-Stillinger algorithm [6]. In this algorithm 
the spheres are compressed uniformly by increasing their 
diameter, while event-driven molecular dynamics is per- 
formed at the same time. We carefully compress the sys- 
tem up to a reduced pressure p = (5Pj ' p = (10 10 -j- 10 12 ) 
in order to easily identify the rattlers. The results ob- 
tained show that already for r = 4, and concentration 
ratio 77 < T] t h, the small component is not blocked, and 
the fraction of rattlers among the small particles is equal 
to 1, see Fig. 2. When this does not happen, the to- 
tal volume fraction of rattlers is always quite small and 
rarely exceeds 5% of the total. For this reason we believe 
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FIG. 2: Fraction of small rattlers fs = N ra ttler,B / Wb. (In- 
set) Volume fraction occupied by the rattlers, defined, follow- 
ing [3], as rjrattier = (/a/»a V 3 (-Da) + f b PbVz(D B )) / , where 
/ M = N ra ttier,y,l 'Nfj, is the fraction of rattlers for species /i. 



that the contribution of rattlers is negligible in this case. 



Additionally, we want to stress that in our theory rattlers 
are not removed from the packing: we assume that their 
cage is vanishingly small so that they are equivalent to 
the jammed particles. This is not a bad approximation 
since, in the region of validity of the theory, the cages in 
which the rattlers are confined are very small. 
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